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ABSTRACT 

We combine deep, wide-field near-IR and optical imaging to demonstrate a reddening- 
independent quasar selection technique based on identifying outliers in the (g — z) 
/ (z — H) colour diagram. In three fields covering a total of « 0.7 deg^ to a depth 
of mH'^18, we identified 68 quasar candidates. Follow-up spectroscopy for 32 objects 
from this candidate list confirmed 22 quasars (0.86<z<2.66), five with significant IR 
excesses. 2 of 8 quasars from a subsample with U band observations do not exhibit 
UVX colours. From these preliminary results, we suggest that this combined optical 
and near-IR selection technique has a high selection efhcicncy (> 65% success rate), 
a high surface density of candidates, and is relatively independent of reddening. We 
discuss the implications for star/galaxy separation for IR base surveys for quasars. We 
provide the coordinate list and follow-up spectroscopy for the sample of 22 confirmed 
quasars. 

Key words: quasars: general - galaxies: active. 



1 INTRODUCTION 

Quasar spectral energy distributions are diverse and com- 
paring surveys at different wavelengths, including the in- 
frared, is important for characterizing the quasar popula- 
tion. In addition, compared to optical surveys, infrared se- 
lection is less effected by dust extinction and reddening (due 
to both dust within host galaxies or along the line of sight). 
This has implications not only for characterizing the influ- 
ence of dust obscuration on the observed population, but 
also for understanding biases in lensing (Kochanek 1996) 
and Damped Lyman-a (DLA) absorption system studies 
based on optical selected quasar samples (Pei and Fall 1995, 
Fall, Pei, McMahon 1989). 

Infra-red observations were first used in quasars selec- 
tion by Braccesi, Lynds and Sandage (1968) where they 
showed how Infrared excess could be used in conjunction 
with UV excess, derived using U— B colours, to distinguish 
quasars from galactic foreground stars. These early obser- 
vations used red sensitive photographic emulsions which 
were sensitive to wavelengths beyond that of the human eye 
but which we would now call the the optical I band(7000- 
9000 A). The underlying physical principle behind the tech- 
nique was the use of a wide wavelength range, 3500-8000A, 



large enough to distinguish the black body dominated stel- 
lar spectra from the non-thermal powerlaw dominated spec- 
tral energy distribution of quasars. The use of IR data 
alone to discover quasars was demonstrated by Beichman et 
al. (1998) who discovered a z=0.147 quasar using JHK data 
from the 2 Micron All Sky Survey (2MASS). However, due to 
current technological limitations, obtaining deep multi band 
IR observations over a sufficiently wide field of view to allow 
the construction of a sample of satisfactory size at cosmolog- 
ically interesting redshifts (i.e. z=l-2) is prohibitively expen- 
sive in terms of telescope time. The availability of large in- 
frared mosaic cameras makes possible wide-field IR surveys 
and the selection of large quasar samples at IR wavelengths. 

The advantage of the IR approach compared to optical 
is the reduced infiuence of both dust extinction and redden- 
ing. The removal of sources from a fiux limited sample by 
dust extinction will be more severe at optical wavelengths 
than in the infrared. In addition, as discussed below, the red- 
dening vector does not push quasars into the stellar locus. 
For example, Francis, Whiting, and Webster 2000 show that 
red quasars are indistinguishable from the stellar locus in 
optical multi-colour surveys, but often can be distinguished 
with the addition of near-IR data. Similarly, Barkhouse and 
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Hall 2001 recommend the combination of optical and near- 
IR as most effective (verse optical or near-IR separately). 

A complementary approach is to use radio or hard X- 
rays as a method to identify AGN in a manner that is not 
strongly biased by dust. In a recent program, Ellison et 
al. (2001) report the results of an survey for DLAs within 
a sample of radio selected quasars. Ellison et al. (2001) 
used a sample of 878 flat spectrum radio source with 
S2.7GHz>250mJy. However it is not obvious that these sur- 
veys can discover enough quasars with z>2, which are lumi- 
nous enough in the optical that searches for damped Lyman- 
a can be carried out, since typically ~1% of quasars have 
radio emission at the 250mJy flux level reached by Elli- 
son et al. (2001). Hooper et al. (1995) observed 3 out of 256 
quasars (1.2%) with mB<18.5 and SgGHz>250mJy , rising to 
14 (5.5%) with S8GHz> 25mJy). Using a deeper radio survey 
is possible but at the mjy level, radio surveys are dominated 
by faint galaxies and the follow-up of such surveys could be 
prohibitive. 

Our approach is similar to the KX technique of Warren, 
Hewett, and Foltz (2000) which proposes selecting quasar 
candidates by their K-excess in a V— J/J— K colour-colour 
diagram (see Groom, Warren and Glazebrook (2001) for an 
application of the method) . As an alternative to V JK based 
candidate selection we demonstrate a technique based on 
gzH observations. There are two reasons for our alternate 
band selection. 

1) Only one near IR band is required. The use of the z 
band, readily obtained with large format GGD based cam- 
eras, leads to an increase in the practicality of surveying 
a large area. The relatively small fleld of view of IR array 
cameras is imposed only on H band observations. 

2) We demonstrate that after accounting for sky bright- 
ness, K-correction and potential differential dust extinction 
effects between the H and K bands, observations in either 
band are equivalent with respect to practical observational 
requirements. The current availability of wide area IR imag- 
ing observations in the H band with GIRSI (section |^ make 
sample deflnition possible over a wide fleld of view (in excess 
of lOdeg^). 

Based on information obtained from the ESO ISSAG 
and SOFI exposure time 



calculator (http://www.eso.org/observing/etc/) we flnd the 
H and K sky brightnesses to be H=14.4 and K=13.0. The 
quasar (H— K) K-correction is estimate to be —0.6 (assuming 
a simple power law quasar model q=— 0.5). Under these as- 
sumptions quasars are 0.8mag brighter relative to the sky in 
H than in K. Strong reddening (Av(rcst)>2.0) in the quasar 
rest frame or along the line of site is required before this sit- 
uation is reversed due to the shallow slope of the reddening 
function at longer wavelengths under most parameterization 
of the reddening law (flgure |^) . 

Although several papers have used the colours of previ- 
ously known quasars to discuss quasar selection using com- 
bined optical and near-IR imaging, we demonstrate this 
technique by carrying out a survey for previously unknown 
quasars. We present follow-up spectroscopic observations 
that conflrm the technique, and suggest high selection effi- 
ciency, high surface density, and relative independence from 
reddening. 

In the following section ^ we discuss the infrared and 
optical imaging data, section (Bl) describes the selection of 
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Figure 1. a) The (H— K) colour for a power law model quasar 
spectrum (q:=— 0.5) is compared to that of a power law with red- 
dening at the level of Av{rost)=l-0- The (H— K) sky colour is 
indicated by the solid line. 

b) Absorption as a function of redshift for Ay =1.0 in the rest 

frame is shown in the gzH and K bands. 

The galactic reddening law of Seaton 1979 is used. 



quasar candidates using colour criteria, section (^ is con- 
cerned with the implications from detection of the host 
galaxy. The spectroscopic observations are described in sec- 
tion (p|), and the sample of confirmed quasars is presented 
in section 

Unless stated otherwise we use conventional Vega mag- 
nitudes and Ho =50 km s^^ Mpc~^, qo=0.5 throughout this 
work. 



2 OPTICAL AND NEAR-IR OBSERVATIONS 

We use optical and near-IR imaging data for three fields 
covering a total of ~ 0.7deg^ to mH~18 (see Table^. These 
fields were chosen due to the availability of both optical and 
near-IR imaging from a number of survey programs. 

The optical imaging was obtained as part of the 
Isaac New ton Telescope Wide Angle Survey (M cMahon et 
al. (2001); http://www.ast.cam.ac.uk/~wfcsur/). This sur- 



vey is being carried out with the prime focus Wide Field 
Gamera (WFG; Ives, TuUoch, and Ghurchill 1996) at the 
2.5m Isaac Newton telescope (INT) on La Palma. The WFG 
consists of a closely-packed mosaic of 4 thinned EEV42 
2kx4k GGDs with a pixel size of 13.5fim corresponding to 
0.33" / pixel and effective field of view of 0.25deg^. Observa- 
tions are taken in 5 wavebands (ugriz) with single exposures 
of 600sec over an area ~100deg^ to nominal 5a limiting mag- 
nitudes of 23, 25, 24, 23 and 22 respectively. The GGD mo- 
saic data is pipeline processed and calibrated in Gambridge 
(Irwin and Lewis 2000). We use a preliminary photometric 
calibration that is accurate to ±0.1 mag. 
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Figure 2. An optical and near-IR colour diagram showing quasar 
model colours (described in the text). The colours of Galactic 
stars are computed from the atlas of Bruzual, Persson, Gunn 
and Stryker included in the IRAF/STSDAS package SYNPHOT. 
The reddening vector does not drive the quasars into the stellar 
locus unlike purely optical colour-colour plots (see figure |ll[) . The 
unfilled circles indicate the colours of known quasars in SDSS 
(Richards et al. 2001), for which we obtained H magnitudes by 
looking for matches (within 2") in the 2MASS second incremental 
release point source catalogue. 



The near-IR imaging was obtained as part of two sur- 
vey projects that use CIRSI (the Cambridge Infrared Survey 
Instrument; Beckett et al. 1996, Mackay et al. 2000). CIRSI 
is a JH-band mosaic imager consisting of 4 Rockwell Ikxlk 
HgCdTe detectors, providing an instantaneous field of view 
of 15.6' X 15.6' on the INT (0.45" / pixel) and 6.6' x 6.6' 
on the du Pont 2.5m (0.19" / pixel) at the Las Campanas 
Observatory. The 1204—0736 field observations were taken 
with the du Pont 2.5m as part of the Las Campanas IR Sur- 
vey (Firth et al. 2001, Chen et al. 2001) during the nights 
29,30 December 1999; 1,2 January 2000; and 14-17 Febru- 
ary 2000. The 1636-1-4101 and 0218-0500 field observations 
were taken with the INT on 13 June 2000 and 4,6,8 Novem- 
ber 1998 respectively as part of the Infra Red CIRSI-INT 
survey. Total exposure times were 2000s for the 1204—0736 
field (up to 4300s in some regions), 2400s for the 1636-1-4101 
field, and 2480s for the 0218-0500 field. The typical obser- 
vation sequence consisted of a 7-point dither with three 30 
second exposures at each dither position, and four telescope 
pointings to fill in the gaps between the chips. 

The CIRSI near-IR data were reduced and mosaiced 
using an automated pipeline (Sabbey et al. 2001), and ob- 
ject catalogs were produced using SExtractor (Berlin and 
Arnouts 1996). The photometric zeropoint was established 
by comparison to 2MASS (see Figure With the precise 
astrometry available in both optical and infrared catalogs 
(0.3" rms with respect to the APM catalog), we were able to 
merge the catalogs by looking for positional matches within 
1". The merged catalog was limited to point sources to re- 
move galaxies from the quasar candidate list. For the spec- 
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Figure 3. The magnitude difference between CIRSI and 2MASS 
catalogues are shown for the 1204—0736 field. The rising tail at 
bright magnitudes is due to saturation in the CIRSI data. The 
median magnitude offset is 0.07m. After correcting for this me- 
dian offset, the rms scatter between the two samples is 0.12m. 
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Figure 4. The observed optical and ncar-IR, colour diagram is 
shown for the 0218-0500 field (Table m). Identified quasars are 
indicated by filled circles. The selection boundary (dot-dashed 
line), chosen based on model quasar colours, is (z — H) > 0.59 X 
(S - 2) - 0.06. 



troscopic observations undertaken to date, the star/galaxy 
seperation was perfromed by requireing a stellar classifica- 
tion for candidate objects in gr and i band imaging data. 
Section |^ discusses the implication of this restriction to the 
quasar candidate list and suggests that the criterion should 
be relaxed to a g band stellar classification alone to prevent 
lower redshift quasars, with resolved host galaxies, being re- 
moved from the candidate list. 
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Figure 5. The follow-up spectroscopy for the quasar sample is 
shown. The first six spectra in the plot were taken with a long- 
slit spectrograph on the du Pont telescope. The remaining spectra 
were taken with the Hydra multi-object spectrograph on WIYN. 
The flux calibrated spectrum of object 6 shows an IRX quasar 
while the fiber spectra of objects 14, 15, 16 and 19 are consistent 
with the IRX photometry. 



QSO Redshift Distribution 




12 3 

RedsPiift z 

Figure 6. The redshift distribution for the sample of 22 con- 
firmed quasars is shown. The observed redshifts range from 
2 = 0.858 to 2 = 2.660 with a median redshift z = 1.211. 



3 COLOUR SELECTION OF Z<3 QUASAR 
CANDIDATES 

Figure |^ shows model quasar colours computed assuming 
an underlying quasar spectrum based on a power law with 
spectral index a=— 0.50 {Sv oc v") and with an emission 
line spectrum based on the Francis et al. (1991) composite 
spectrum. The influence of a range of values for a is in- 
dicated in the plot. The reddening vector is derived over 
a range of redshifts using the galactic reddening model of 
Seaton (1979). Synthetic photometry calculations are per- 
formed with the IRAF/STSDAS package SYNPHOT. The 
stellar locus is calculated syntheticaly using the stellar atlas 
of Gunn-Stryker (1983) using the extended spectral cover- 
age spectra from the Bruzual-Persson-Gunn-Stryker spec- 
trophotometric atlas found in SYNPHOT. Representative 
filter transmission curves for the observing system are ap- 
plied to the input spectra to predict the colours of objects. 

The optical and near-IR colour diagram (flgure|^ shows 
that quasar candidates are expected to lie separated from 
the stellar locus in colour space and that reddening due to 
dust, either in the host galaxy or along the line of sight to the 
quasar, will not drive the quasars into the stellar locus (in 
contrast to purely optical colour-colour plots, see figure |ll[ ). 
Reddening of foreground stars moves stars parallel to the 
stellar locus preventing a broadened stellar locus overlapping 
the quasar region. 

The mainly optically selected quasar sample of 
Richards et al. (2001) with SDSS photometry and matched 
with the 2MASS second incremental release is over- 
layed for comparison. There is evidence that the 
colours of the three quasars hidden with in the stel- 
lar locus in this plot are strongly effected by vari- 
ability (NED identiflers NGC0450:0126-0019 (z=1.76), 
UM357:0140-0050 (z=0.33), PKS1215-002 z=(0.42)). One 
object, UM357:0140-0050, is a known Optical Violent Vari- 
able (OVV) and all three objects are located on the quasar 
power law spectrum locus in g— z v J— K colour space. This 
latter space is insensitive to variability since for 2MASS 
the IR data is recorded simultaneously through the use of 
dichroics and for the SDSS bands the lag between obser- 
vation is insignificant (~minutes) when compared to the 
variability timescale for broad band optical observations of 
quasars. 

To identify quasar candidated the stellar locus is ap- 
proximated by a two component linear fit of the form. 



{z — H) > m X 



z) + c + offset 



(1) 



The offset value being chosen to minimize the number of 
stray stars in the sample while rejecting a minimum of 
quasar candidates. 

The region above these boundaries (lower g— z for a give 
value of z— H) defines the candidate selection region. All stel- 
lar objects in this region are visually inspected to check for 
spurious photometry due for example to stellar diffraction 
spikes or cosmic rays in the single exposure INT WFS data. 

At redshift z>3, the Lyman-a forest passes into the g 
band and absorption in the Inter Galactic Medium (IGM) 
moves quasar colours into the stellar locus. Few z>3 quasars 
would therefore be expected in a gzH selected sample. 
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3.1 Prospects for extending the technique to 
z>3.0 

As part of the Isaac Newton Telescope Wide Angle Survey 
(Walton, Lennon, Irwin and McMahon 2001) optical obser- 
vations, over the full range of wavelengths ugri and z, are 
available for two of the fields examined in this paper. This 
data set has been successfully used for the identification of 
high redshift (z~5) quasars via optical colour selection tech- 
niques (Sharp, McMahon, Irwin and Hodgkin 2001). How- 
ever, the gzH technique can not be used at higher redshift 
(z>3, after the onset of IGM absorption in the g band) by 
replacement of the g band alone with the r or i bands as 
quasars merge with the stellar locus (as shown in the colour 
diagrams of figure ^). Figure ^ shows, as discussed by Bark- 
house and Hall (2001), it is possible to use the RJK diagram. 
Pushing the technique into the i band and using the iJK dia- 
gram allows the selection of targets out to redshift 5 (before 
IGM absorption confuses the selection with the low mass 
region of the stellar locus) but blending with the main se- 
quence is greater than that in the rJK diagram. Figure ^ 
demonstrates a range of possible selection plots. 

The lower surface density of 3<z<4 quasars, in compar- 
ison to the high selection efficiency (65%) gzH z<3 quasars, 
implies that a rzH selected quasar sample could be effec- 
tively integrated into a joint IR quasar selection campaign, 
using wide field of view multi-object spectrograph observa- 
tions, capitalizing on the high multiplex capabilities of in- 
struments such as HYDRA. 

4 POTENTIAL EFFECTS OF THE HOST 
GALAXY 

There are two potential problems associate with the detec- 
tion of the quasar host galaxy. Both are illustrated in figure 

i 

1) The magnitude limit of the gzH selection technique 
is essentially set by the depth of the H band data due to 
the relative easy with which deep g and z observations can 
be achieved. However the depth of the INT WFS imaging 
data is sufficient that potentially the host galaxy of quasars 
may be detected and spatially resolved in one or more of the 
optical passbands. Star/galaxy separation is used to remove 
many galaxies from the sample and this step may also reject 
lower luminosity quasars with bright host galaxies. 

2) The addition of a galaxy component to the observed 
flux may alter the objects colour to the extent that it would 
no longer fall within the colour selection boundary. 

The absence of any quasar with z<0.8 in our current 
sample is indicative of the first problem. Using the lumi- 
nosity function of Boyle et al. (2000) , and within the uncer- 
tainties of the conversion from mn to Mb, we expect ~2.6 
quEisars per deg^ to mH<18.0 in the range 0.0<z<0.8 and 
Mb<— 23. However with only ~47% of candidates observed 
at this time and a nominal survey area of ~0.7deg^ this be- 
comes ~0.9 quasars that may be missing from the sample, 
not inconsistent with the small sample size. 

Above z=0.8 the optical g band, with which star/galaxy 
separation is performed, samples the rest frame UV contin- 
uum (longword of the 4000A break) where the host galaxy 
presence in imaging data will be at it's weakest. 
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Figure 8. The magnitude evolution tracks are shown for an early 
type galaxy and the model quasar spectrum. Both objects are 
normalized to M*=MBj=-19.5 (Madgwick et al. 2001) for the 
comparison. Figure ^ demonstrates the observed range in quasar 
and host galaxy magnitudes. The galaxy template used is taken 
from Coleman, Wu and Weedman 1980 and no evolution in the 
spectrum is invoked. The nominal survey 5cr limits are marked by 
the dashed horizontal line. While the galaxy contribution to the 
g band is weak it may dominate the flux in the IR H band. The 
shaded area indicates a 2 magnitude wide zone below the PSF 
limited quasar magnitude. If the host galaxy aperture magnitude 
lines within this zone an object may be classified as an extended 
source. 



Morphological classification of objects within the INT 
WFS is based on an analysis of the curve of growth of an 
objects fiux in apertures of radius 1/2, 1 ,y/2, 2, 2y/2 times 
the seeing for the observation. A detailed discussion of the 
classification algorithm is beyond the scope of this paper. 
However investigation of the aperture flux ratio and classi- 
fication of simulated and observed early type galaxies and 
PSF limited point source profiles show that quasar cores will 
be classified as stellar in the INT WFS g band in the pres- 
ence of an extended host galaxy component provided the 
host is 2 magnitudes fainter than the quasar core (assum- 
ing a De Vaucouleur light profile for the host galaxy with 
half-light radius 0.5" and axis ratio 0.8). This comparison is 
based on the seeing radius aperture magnitude as oppose to 
a total magnitude. A stellar classification will be obtained 
for the quasar core of a composite object for a lower sepa- 
ration in total magnitude, the value depending strongly on 
the definition of total magnitude. 

Figure |l^ addresses the second effect. A reproduction 
of figure ^ is shown with the colour evolution track (0<z< 
3, 0.1 steps in z) for an non evolving Early type galaxy 
overlayed. A colour evolution track assuming a combination 
of the model quasar and such a host galaxy, with the galaxy 
normalized to 4 magnitudes fainter than the quasar in the 
rest frame B band, is indicated by the solid line. 

While it is clear that any host galaxy contribution to 
the quasar colour has a significant effect, it acts to increase 
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Figure 7. A range of potential optical-IR quasar colour selection diagrams are considered. A sample of quasars is constructed from the 
SDSS quasar sample of Richards et al. (2001) with IR data obtained from a 2" radius match to the 2MASS second incremental release 
data. A stellar locus is computed from the spectral atlas of Gunn and Stryker (1983). Reddening vectors are shown for Av(rest)=1.0 
under the Seaton (1979) galactic reddening law. The quasar locus with redshift is marked for a range of redshifts. 



© 0000 RAS, MNRAS 000, 000-000 



Reddening- Independent Quasar Selection 7 




Figure 9. A comparison of Mb for a sample of quasars, taken 
from tlie LBQS sample and assuming elliptical host galaxies, is 
obtained from Hooper, Impey and Foltz (1997). Conversion from 
Mr to Mb is performed using the quasar model described in the 
text and the galaxy template of Coleman, Wu and Weedman 1980. 
Values for Mr where derived from profile fitting and represent 
total magnitudes. 



Figure 10. The gzh diagram from figure |2| is shown overlayed 
with the colour evolution track for a non evolving early type 
galaxy (0<z<3, 0.1 z step). The colour track for a combination of 
quasar and elliptical galaxy spectrum is shown with the elliptical 
component normalized the 4 magnitudes below the quasar in the 
rest frame B band. The galaxy component dominates the H band 
flux. 



the detectability of the quasar in the sense that it improves 
separation from the stellar locus. 

From Figure ^, the nominal H band magnitude limit 
of our survey of mn^lS.O corresponds to an absolute B 
band magnitude of —24.5 and —25.5 at z=l and 2 respec- 
tively. Therefore if the quasar host galaxy is uncorrelated 
with the quasar luminosity and is assumed to correspond 
to an unevolved L* ellipical galaxy (Mb =—19.5, Madgwick 
et al. 2001) the quasar host will be 3-5 magnitudes fainter 
than the quasar nucleus. However, if the quasar host galaxy 
is correlated with the quasar luminosity such as one would 
expect from the observed correlations in the local Universe 
between dormant quiescent Black hole masses and galaxy 
bulge luminosities (Magorrian et al. 1998) the quasar host 
galaxy may contribute to the H band light. Taking the ob- 
served correlation between rest frame absolute magnitudes 
of z~0.5 quasars and quasar hosts observed by Hooper, Im- 
pey and Foltz (1997) where AMb=2 (see figure|^ the quasar 
host may contribute ~50% of the hght in H in the redshift 
range 1 to 2 (see figure 



blue end was set by a long pass filter). The dispersion was 
2.64A/pix and the resolution was 5.7A. Arc lamp (Cu-Ar) 
calibration frames were taken to set the wavelength scale. 
With the 1 degree field diameter and 97 red fibers (2" diame- 
ter) available, we were able to obtain spectra for both quasar 
candidates and a sample of random targets spread through- 
out the colour diagram. The minimum exposure time was 
1.7 hours, although many objects were observed twice in 
overlapping Hydra fields resulting in 3.4 hour exposures. 

Long-slit spectroscopic observations were taken during 
the five nights from 22 to 26 February 2001 at the du Pont 
2.5m telescope at Las Campanas Observatory. The weather 
conditions and seeing were both good. The Modular Spec- 
trograph was used with the SITe#2 detector and the 150mm 
grating to obtain spectra with 4A/pixel over the wavelength 
range 4300-8500A. We obtained spectra for 38 objects (in- 
cluding separate projects), with exposure times ranging from 
10 minutes to 1 hour. Wavelength calibration was performed 
with arc spectra. 



5 SPECTROSCOPIC OBSERVATIONS 

Multi-object spectroscopy was obtained using the Hydra 
multi-fiber spectrograph (Harden and Armandroff 1995) on 
the 3.5m WIYN telescope at the Kitt Peak Observatory. The 
1637+4101 and 0218-0500 fields were observed during 2 and 
3 September 2000 (as part of an observing program contain- 
ing other projects) in clear sky conditions with a seeing of 
~ 0.8". We used the Bench Spectrograph Camera and a 316 
lines/mm grating with a central wavelength of 6027A. The 
total wavelength coverage was from 4240A to 8728A (the 



6 CONFIRMED QUASARS 

The confirmed quasars are listed in Table ^ Object 16 has a 
possible broad emission-line at 6500A, but the spectrum is 
low signal-to-noise and requires independent confirmation. 
Based on Veron-Cetty and Veron (2000), 2 of the 22 con- 
firmed quasars were previously known. Both known quasars 
where discovered as part of the slitless spectroscopy survey 
of Crampton et al. (1988). Three additional objects are listed 
in the NED data base as candidate quasars from Crampton, 
Cowley, Hartwick and Ko (1992) but without prior spectro- 
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Figure 11. Multi wavelength data is available in ugrizJH for a 
subset of the current sample. Detailed analysis of the multi band 
properties of the full quasar sample is underway. Here we show 
a comparison of a sample of previously known quasars identified 
in SDSS (Richards et al. 2001) with the new quasars identified by 
gzH selection for which ugr data is available. The horizontal line 
at u— g=— 0.3 represents a UVX selection boundary analogous to 
that used in the 2df quasar survey (2QZ Boyle et al. 2000). The 
stellar locus, computed from Gunn and Stryker (1983), is shown 
along with the observational data from one field (~0.25deg2) from 
the INT WFS. The locus of quasar colour with redshift is indi- 
cated by the solid line. Reddening vectors for dust models based 
on Galactic (Gal-Seaton 1979) and Small Magelanic Cloud (SMC- 
Prevot et al. 1984) are shown over a rang e of redshifts. The 2200A 
feature in the galactic law passes through the g band over the 
range l<z<2. 



scopic confirmation and redshifts. The redshift distribution 
is shown in figure ^. The median redshift is z=1.211. 

The current spectroscopic sample size hampers a full 
comparison of the quasar space density with other works. 
The sample is mn limited making comparisons with optical 
samples complex due to uncertainties in the K-correction 
required to compute optical absolute magnitude indicators. 
The sample is not complete to an optical magnitude limit in 
mg or mz due to incomplete follow up of the fainter optical 
candidates. 

U band observations are available for a subsample of 8 
quasars. Two quasars, 17 and 22, do not exhibit strong UV 
emission and could be missing from a classical UVX selected 
sample. Defining a IRX criterion of z-H>1.5, neither object 
is IRX (z-H=1.27 and 1.37 respectively). However, five IRX 
quasars are identified in the sample, quasars 6, 14, 15, 16 
and 19. Objects 15 and 16 show low UVX (-0.35 and -0.30) 
while quasar 19 has a more typical value of -0.6. No u band 
observations are available for objects 6 and 14 at this time. 

Figure |l^ compares the histogram of quasar colours 
for the 22 quasar CIRSI sample and 89 redshift z<2.6 
quasars from the Richards et al. (2001) sample matched to 
the 2MASS point source catalogue to provide H band mag- 
nitudes. The Richards et al. (2001) sample is normalized to 
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Figure 12. A comparison of the histograms of quasars colours 
for the CIRSI sample and that of Richards et al. (2001), from 
SDSS and matched to the 2MASS point source catalogue, shows 
broad similarity between the distributions. The SDSS sample (89 
quasars with z<2.6) is scaled to the CIRSI sample size of 22 
quasars. A K-S test between the distributions yields probabilities 
that they of are drawn from the same parent population of : 
K-z=0.82, z-H=0.76. 



the total sample size of the CIRSI observations. Within the 
limited sample size the colour distributions are similar with 
a possible excess of IRX quasars in the CIRSI sample. 

As discussed in section ^, the lack of quasar in the sam- 
ple with redshift z<0.8 is not inconsistent with the luminos- 
ity function predictions within the constraints of the current 
survey area (0.7deg^) and ~47% candidate follow up. The 
low number of quasars at z>2 is also not inconsistent. 
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Figure 13. The redshift z=0 quasar and elliptical galaxy spectra 
used in the discussion are shown against the rest frame filter pass 
bands. Both SEDS are normalized to 1 at A4400A , the efli'ective 
wavelength of the B bands. The optical ugriz and IR JHK filter 
sets are shown for comparison. Quasar - S,^ oc u°' (o=— 0.50) with 
an emission line spectrum based on the Francis et al. (1991) com- 
posite spectrum , Elliptical galaxy - Coleman, Wu and Weedman 
(1980). 



7 CONCLUSIONS 

We report the identification of a gzH selected sample of z<3 
quasars. The g — z / z — H colour diagram is used with high 
selection efficiency (~65%) and reduced sensitivity to red- 
dening and extinction. The initial sample has 22 confirmed 
quasars, but the full INT CIRSI survey will contain several 
hundred. To date two objects within a subsample of the iden- 
tified quasars, for which UV observations are available, do 
not possess a significant UV excess and would not have been 
recovered by a UVX based selection. The full CIRSI-INT 
gzH quasar sample will cover an area of in excess of lOdeg'^ 
including observations across the full range of bands ugrizj 
and H. Contrasting the sample with those compiled through 
a range of selection techniques, such as UVX colour selec- 
tion and radio identifications from overlapping surveys such 
as the VLA-FIRST survey, will provide an excellent tool for 
identifying biasing within quasar selection methodologies. 
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Table 1. Survey fields 



Name R.A. (J2000) Dec Area Depth (Vega) 











deg^ 


9 


2 


H 


1204- 


-0736 


12:04:50 


-07:36:00 


0.18 


24 


21 


18 


1636- 


h4101 


16:36:50 


+41:01:50 


0.35 


24 


21 


18 


0218- 


-0500 


02:18:00 


-05:00:00 


0.17 


24 


21 


17.5 



Table 2. The quasar Sample. 



Name 



R.A. (J2000) Dec 



magnitude (Vega) 



Comments 











nig 


rriz 


van 










CIRSIl 


12 


03:48.61 


-07:19:00.7 


19.98 


18.76 


17.77 


0.858 


1.22 


0.99 




CIRSI2 


12 


04:21.62 


-07:22:01.0 


19.46 


18.41 


17.78 


1.189 


1.05 


0.63 




CIRSI3 


12 


04:09.20 


-07:24:22.8 


19.89 


18.96 


18.12 


1.625 


0.93 


0.84 




CIRSI4 


12 


04:08.18 


-07:21:45.6 


20.13 


19.41 


18.26 


1.097 


0.72 


1.15 




CIRSI5 


12 


05:34.59 


-07:41:21.4 


19.08 


17.97 


17.19 


2.660 


1.11 


0.78 




CIRSI6 


12 


05:07.68 


-07:45:26.6 


21.01 


19.40 


17.89 


1.420 


1.61 


1.51 


IRX, No UV data 


CCS88 163351.6-1-410628 


16 


35:31.05 


-1-41:00:27.2 


19.14 


18.25 


17.21 


1.151 


0.89 


1.04 


Crampton et al. (1988) 


CIRSI8 


16 


37:10.03 


+40:56:42.9 


20.21 


19.22 


18.06 


1.434 


0.99 


1.16 


Candidate Crampton et al. (1992) 


CIRSI9 


16 


37:00.64 


+41:05:55.2 


19.69 


18.78 


18.10 


2.061 


0.91 


0.68 


Candidate Crampton et al. (1992) 


CIRSIIO 


16 


36:47.16 


+41:03:35.0 


20.54 


19.37 


18.29 


1.077 


1.17 


1.08 




CIRSIl 1 


16 


37:14.80 


+41:12:32.6 


20.16 


19.18 


18.50 


1.642 


0.98 


0.68 


Candidate Crampton et al. (1992) 


CIRSI12 


16 


35:30.51 


+41:10:41.6 


20.86 


19.36 


18.04 


1.211 


1.50 


1.32 




CCS88 163447.3-1-405448 


16 


36:27.11 


+40:48:48.9 


19.72 


18.89 


17.90 


0.904 


0.83 


0.99 


Crampton et al. (1988) 


CIRSI14 


16 


37:34.03 


+41:16:09.0 


22.42 


20.05 


18.27 


1.384 


2.37 


1.78 


IRX, No UV data 


CIRSI15 


02 


17:56.51 


-05:06:51.7 


20.72 


19.79 


18.01 


1.087 


0.93 


1.78 


IRX, low UVX 


CIRSI16 


02 


18:34.44 


-05:13:56.9 


20.95 


19.04 


17.32 


1.351 


1.91 


1.72 


IRX, low UVX 


CIRSI17 


02 


17:11.99 


-04:46:19.8 


18.92 


18.00 


16.73 


1.102 


0.92 


1.27 


Non UVX 


CIRSI18 


02 


18:08.56 


-05:12:23.9 


20.22 


19.43 


18.18 


1.038 


0.79 


1.25 




CIRSI19 


02 


18:17.42 


-04:51:12.5 


19.41 


18.19 


16.68 


1.085 


1.22 


1.51 


IRX 


CIRSI20 


02 


18:13.94 


-04:52:41.1 


19.80 


18.67 


17.43 


1.443 


1.13 


1.24 




CIRSI21 


02 


18:30.57 


-04:56:22.7 


17.60 


16.61 


15.61 


1.401 


0.99 


1.00 




CIRSI22 


02 


17:21.66 


-05:06:28.9 


19.97 


18.68 


17.31 


0.983 


1.29 


1.37 


Non UVX 
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